Collective Sliding States for Colloidal Molecular Crystals 
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We study the driving of colloidal molecular crystals over periodic substrates such as those created 
with optical traps. The n-merization that occurs in the colloidal molecular crystal states produces 
a remarkably rich variety of distinct dynamical behaviors, including polarization effects within the 
pinned phase and the formation of both ordered and disordered sliding phases. Using computer 
simulations, we map the dynamic phase diagrams as a function of substrate strength for dimers 
and trimers on a triangular substrate, and correlate features on the phase diagram with transport 
signatures. 

PACS numbers: 82.70.Dd 



A wide range of condensed matter systems can be mod- 
eled as a collection of interacting particles on a periodic 
substrate under equilibrium and nonequilibrium condi- 
tions, where the number of particles may be commensu- 
rate or incommensurate with the substrate. Examples of 
such systems include molecules and atoms on corrugated 
surfaces [H, [3], models of sliding interfaces [3], vortices 
interactingwith artificial pinning arrays in superconduc- 
tors 0, 0, Q , and vortices in Bose- Einstein condensates 
interacting with optical trap arrays f^. Recently, there 
has been growing interest in studying colloidal particles 
interacting with one-dimensional [8| or two-dimensional 
periodic or quasiperiodic substrates such as those gener- 
ated with optical traps. Sufficiently strong traps 
can capture multiple colloids, and at commensurate fill- 
ings where the number of colloids is equal to an integer 
multiple of the number of traps, the colloids in each trap 
undergo an effective n-merization which causes them to 
act like rigid objects such as dimers or trimers. This 
produces an orientational degree of freedom and leads to 
the formation of ordered states which have been termed 
colloidal molecular crystals 14, 3, 3 III- Simulations 
and experiments have shown that these systems can un- 
dergo interesting ordered to orientationally ordered to 
disordered transitions as a function of substrate strength 
14, . Additional theoretical and numerical studies 
demonstrated that colloidal molecular crystal systems 
can exhibit a rich variety of equilibrium states which 
have ferromagnetic, antiferromagnetic, and other types 
of spinlike symmetries 1^ 17 1. 

Although individual colloids have been driven over pe- 
riodic substrates flll,fl^, and colloidal assemblies have 
been depinned from random substrates [l3| , the nonequi- 
librium dynamics of colloidal molecular crystals in the 
presence of an additional driving force has not been stud- 
ied previously. The n-merization that occurs in the col- 
loidal molecular crystal system could be expected to pro- 
duce dynamical sliding behaviors that are distinct from 
those of sliding point particles. In this work, we show 
that the depinning force passes through a series of peaks 
as a function of colloid density that are associated with 



the formation of commensurate colloidal molecular crys- 
tal states. At the commensurate densities, the scaling 
of the depinning threshold indicates that the pinning is 
collective, while pronounced changes occur in the depin- 
ning threshold as a function of substrate strength that 
correlate with changes in the symmetry of the colloidal 
molecular crystal states. Changes in the structure of the 
colloidal molecular crystals also produce features in the 
velocity-force curves. We specifically map the dynami- 
cal phase diagram for dimer and trimer states at filling 
fractions of 2 and 3 colloids per trap. In addition to pro- 
viding an understanding of the general phenomenon of 
sliding phases on periodic substrates, our results could 
also be of importance for the development of externally 
driven dynamical assembly techniques for colloids and 
other particulate matter systems. 

We consider a two-dimensional system with periodic 
boundary conditions in the x and y directions containing 
Nc colloidal particles. The dynamics of a single colloid i 
at position is governed by the overdamped equation 
of motion I™ 



dt 



(1) 



where we set 77 = 1. The colloid-colloid interaction force 
is F^'^ — —Y^^^i'^V{Rij) where the potential has a 
Yukawa form, V{Rij) — {Eq/ Rij) exp{~KRij) and where 
i?y = |Ri - Rjl, Eq = Z*'^ /(ATreeoao), e is the solvent 
dielectric constant, Z* is the effective charge, and 1/k is 
the screening length. Lengths are measured in units of 
ao, assumed to be on the order of a micron, forces are 
measured in units of Fq = Eg/ao, and time is measured 
in units of t = ti/Eq. The substrate force Fg arises from 
a triangular substrate with 



sm 



2TTb, 

ao 



[cos(0,)x-sin(0,)y], (2) 



where bi = x cos{6i)--y sm{9i)+ao /2, 9i — 7r/6, 6*2 — 7r/2, 
and 6*3 = 57r/6. Here A is the relative substrate strength 
and there are 7V« substrate minima. The initial colloidal 
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FIG. 1: (a) The critical depinning force Fc versus filling 
fraction Nc/Ns for a substrate with strength A = 2.5. Peaks 
appear at the commensurate fillings of Nc/Ns — 2, 3 and 
4. (b) V, the average colloidal velocity per particle, vs the 
external drive Fd for Nc/Ns = 1.94 (dotted fine), 2 (thick 
line), and 2.13 (thin line). A clear depinning threshold exists 
which is largest at Nc/Ns = 2. (c) The colloidal positions 
(black dots) and trajectories (black lines) for Nc/Ns — 2.0 in 
the pinned triangular crystal (PC) state at Fd = 0.0 and A = 
0.5. (d) Colloid positions and trajectories for Nc/Ns = 2.0 in 
the moving random (MR) state at A — 1.5 and Fd/Fc — 1.1. 
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FIG. 2: The colloid positions (black dots) and trajectories 
(black lines) for Nc/Ns = 2.0 and (a-d) A = 3.25; (e,f) 
A — 1.5. (a) The pinned herringbone (PHB) state at Fd — 0. 
(b) The PHB at finite Fd/Fc = 0.45 showing the onset of 
dimer polarization in the direction of the drive, (c) The fully 
polarized dimers at Fd/Fc — 0.9 form the pinned ferromag- 
netic (PF) state, (d) The moving ferromagnetic (MF) state 
at Fd/Fc = 1.1 where the motion is strictly one-dimensional, 
(e) The PMF state for A = 1.5 at Fd/Fc = 1.5. (f) Vortex 
positions only in the PMF state from panel (e) showing that 
every other row of dimers is aligned. 



positions arc obtained through simulated annealing. The 
applied driving force = FdX represents the force that 
would be produced by an external electric field [l3| . For 
each drive, we measure the average colloid velocity V = 
{N~^ dRi/dt ■ x). The depinning force Fc is defined 
as the value of at which V = 0.025. 

In Fig. [Tfa) we plot the depinning threshold Fc versus 
the filling factor Nc/Ng over the range 1.5 < Nc/Ns < 4.5 
for a system with a substrate of strength ^ = 2.5. At the 
integer matching fields Nc/Ng — 2, 3, and 4, there are 
clear peaks in the depinning threshold. This is similar 
to the behavior observed for vortices in superconductors 
with periodic pinning, where peaks in the critical cur- 
rent (which is proportional to the depinning force) ap- 
pear when the vortex density is an integer multiple of 
the pinning site density 0, 0] . In Fig. (Ub) we illustrate 
the velocity force curves for Nc/Ns = 1.94, 2.0, and 2.13. 
A clear single depinning threshold appears at Nc/Ns = 2, 
while two-step depinning transitions occur at the nonin- 
teger fillings. For Nc/Ns — 1.94, the initial depinning oc- 
curs due to the motion of monomer defects in the dimer 
lattice followed by the depinning of the remaining dimers, 
while for Nc/Ns — 2.13, the trimer defects in the dimer 
lattice depin first and then the remaining dimers depin 
at a higher drive. 

We measure Fc versus substrate strength A for the 
case Nc/Ns = 2.0 and find that there are three distinct 
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FIG. 3: (a) V vs Fd for Nc/Ns = 2.0 and A = 3.25 showing 
the contirmous PF-MF depinning transition. Inset: Vv&Fd — 
Fc in the same system indicating a power-law scaling with a = 
0.5. (b) V vs Fd for Nc/Ns = 2.0 and A = 1.5. The depinning 
occurs from a PHB state to a moving random (MR) state, 
and at higher drives the colloids organize into a partially- 
ordered moving ferromagnetic (PMF) state. Inset: dV/dFd 
vs Fd for the same system. The sharp jump indicates that the 
depinning transition is discontinuous. The fluctuating region 
corresponds to the MR state and the onset of a regime with 
small fluctuations at _fd > 1 corresponds to the formation of 
the PMF state. 



pinned states and four moving states. In Fig. [5] we plot 
the colloidal configurations at different values of Fd for 
A = 3.25. At Fd = 0, the ground state is the pinned 
herringbone (PHB) structure shown in Fig. [H^a). For 
increasing Fd, the dimers become increasingly polarized 
in the direction of drive, as illustrated in Fig. [D^b) for 
Fd/Fc = 0.45. At sufficiently large Fd, the dimers are 
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completely aligned into the pinned ferromagnetic (PF) 
state shown in Fig.[2l^c) for Fd/Fc = 0.9. The system de- 
pins directly into the moving ferromagnetic (MF) state, 
where the colloids move in one-dimensional channels as 
seen in Fig. ^d) . This same sequence of states also ap- 
pears for higher values of A. The PF-MF transition is 
continuous, as shown in Fig.[3l^a) where we plot V versus 
Fd for A = 3.25. The inset of Fig. [Sja) indicates that 
the velocity follows a power law scaling, V — {Fd — Fc)" 
with a = 1/2, with a turnover at higher Fd to a linear 
form. This exponent is consistent with elastic depinning, 
in which the colloids keep the same neighbors while mov- 
ing. We observe the same velocity scaling for A > 2.5. 

For 0.7 < A < 2.5, the PF state does not form. In- 
stead, the PHB depins discontinuously into a plastically 
flowing state where the colloids can exchange neighbors 
and undergo a transverse diffusive behavior shown in 
Fig. [Hd). We call this the moving random (MR) state. 
In Fig. [3fb) we plot the velocity force curve a,t A = 1.5 
which has a discontinuous depinning transition, as indi- 
cated in the inset of Fig. [SJb) where we show dV/dFd 
versus Fd. There is a clear sharp jump in V at the 
PHB-MR depinning transition, followed by a regime of 
fluctuating dV/dFd which corresponds to the MR state. 
For Fd > 1.0, the fluctuations are diminished when the 
system forms a partially-ordered moving ferromagnetic 
(PMF) state, illustrated in Fig. ^e) and (f). Here the 
moving colloids form a stripelike structure in which every 
other row of dimers is aligned with the drive. Unlike the 
MR state shown in Fig.[ljd), in the PMF state there is no 
transverse diffusion. For A < 0.7 and Fd = 0, the weak 
substrate causes the elastic energy cost of the PHB state 
to be too high, and instead the colloids form a pinned tri- 
angular crystal (PC) state illustrated in Fig. [ijc). The 
PC state is only weakly pinned and the depinning occurs 
as a continuous elastic transition to a moving triangular 
crystal (MC) state. 

In Fig. m we summarize the different states in a dy- 
namic phase diagram for Nc/Ns ~ 2.0. We plot dF^/dA 
versus A in the upper inset of Fig. |3]to show how changes 
in Fc correlate with different phases of the system. The 
PHB-PC transition is marked by a discontinuous jump 
down in Fc when the PC phase forms, as indicated by the 
sharp peak in dF^/dA. A smaller feature in dF^/dA oc- 
curs at A = 2.5 when the depinning changes from the dis- 
continuous PHB-PMF transition to the continuous PF- 
MF transition. The PF phase forms for A > 2.5 and 
the PHB-PF transition shifts to lower Fd with increas- 
ing A. Since the PHB state forms due to the effective 
quadrupole interaction between the dimers [lB|, as A in- 
creases, the colloids forming each dimer are pulled closer 
together, reducing the quadrupole moment and facilitat- 
ing the formation of the PF and MF states. The MR state 
appears in a narrow window between the two elastic de- 
pinning transitions, PHB-MC and PHB-PMF. In the MR 
regime, there is a competition between the quadrupole 




FIG. 4: The dynamic phase diagram for Fd vs A with 
Nc/Ns — 2.0. Open circles: depinning threshold. Dashed 
line: a fit to a power law with a — 2.0. PC, pinned triangular 
crystal; PHB, pinned herringbone; PF, pinned ferromagnetic; 
MC, moving triangular crystal; MR, moving random; MF, 
moving ferromagnetic; and PMF, partially-ordered moving 
ferromagnetic. Upper left inset: dFc/dA vs A for the depin- 
ning curve in the main panel. The sharp peak separates the 
PC and PHB states. Lower right inset: Fc vs A for a single 
isolated particle, showing a linear scaling. 



moment, which prevents the dimers from aligning into 
the PF state, and the tendency of the external drive to 
align the dimers. For 0.7 < A < 1.25, the MR state or- 
ders into a MC state for increasing Fd rather than form- 
ing the PMF state. For small A, Fc exhibits a scaling 
Fc oc A^ , consistent with collective pinning in two di- 
18|. A similar scaling occurs for A > 0.7 as 



mensions 

well. In contrast, for a single colloid moving over the pe- 
riodic substrate, we find the linear behavior Fc cx A, as 
shown in the lower inset of Fig. 21 The collective pinning 
behavior that occurs in the PHB state for large A likely 
arises because the objects that are depinning are dimers 
rather than single particles. 

For Nc/Ns = 3.0, we find a similar set of pinned and 
moving phases as shown in the dynamic phase diagram 
for Fd vs A in Fig.[5]Ja). At this filling, there is only a sin- 
gle pinned phase, the pinned trimer crystal (PTC) illus- 
trated in Figs.in^b) and (c) for A = 0.5 and A — 2.5, and 
so there is no discontinuous transition in Fc vs A such as 
that seen at the PHB-PC transition for Nc/Ns = 2.0. In 
Fig.ElJa), the PTC-MC depinning transition for A < 1.0 
is elastic, while for A > 1.0 the PTC-MR depinning tran- 
sition is followed at higher drives by the organization of 
the colloids into either the moving crystal (MC) state or 
into the moving stripe (MS) state illustrated in Fig.[5l^d). 
In the MS state, the colloids are strictly confined to move 
along one-dimensional rows with no transverse diffusion. 
The oriented trimer structure found in the PC is lost in 
the MS state; however, some orientational ordering of 
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FIG. 5; (a) The dynamic phase diagram for Fd vs A with 
Nc/Ns = 3.0. Open circles: depinning threshold. PTC: 
pinned trimer crystal; MC: moving triangular crystal; MR: 
moving random; MS: moving stripe, (b) Colloid positions in 
PTC at ^ = 0.5 and Fa = 0. (c) Colloid positions in PTC 
at A = 2.5 and Fd = 0. The size of each trimer is reduced 
compared to panel (b). (d) The MS state at ^4 = 2.5 and 
Fd = 1.5. 



the trimers persists in the MS state, producing a zig- 
zag structure. The MR state reaches its maximum ex- 
tent at the transition between the MC and MS phases 
a.t A = 1.6. The scaling of Fc with A is consistent with 
collective depinning for low A, while for A > 3.0, there 
is a rollover to a more linear regime, indicative of single 
particle depinning. In the high-A regime, the depinning 
occurs via the hopping of individual colloids from well to 
well, followed at higher drives by more general motion of 
all the colloids. 

We expect that a similar type of phase diagram will 
occur for the higher order molecular crystals. Addi- 
tional phases are likely to occur at incommensurate fill- 
ings. Thermal fluctuations could produce interesting re- 
sults and new effects since it has been shown that the 
disorde ring transition depends on the substrate strength 

In summary, we find that a remarkably rich variety 
of dynamical sliding states can be realized for colloidal 
molecular crystals under an external drive with varied 
substrate strength. For dimer colloidal molecular crys- 
tals, these states include moving ordered, partially or- 
dered, and moving random phases. The external drive 
can induce transitions within the pinned phase itself such 
as a polarization from a pinned herringbone to a pinned 
ferromagnetic state. The onset of the different phases 
can be identified through features in the transport re- 
sponse and depinning threshold. We map the dynamical 
phase diagrams for dimer and trimer states and find that 
similar features exist for both fillings. Our results should 



be useful for the general understanding of sliding states 
for complex particles on periodic substrates and may also 
have potential applications for creating dynamically in- 
duced self-assembled structures. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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